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Abstract 

We show that it is possible to devise a large class of skew-product dynamical 
systems which have strange nonchaotic attractors (SNAs): the dynamics is 
asymptotically on fractal attractors and the largest Lyapunov exponent is 
nonpositive. Furthermore, we show that quasiperiodic forcing, which has 
been a hallmark of essentially all hitherto known examples of such dynamics 
is not necessary for the creation of SNAs. 
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I. INTRODUCTION 



Since 1984, when Grebogi, Ott, Pelikan, and Yorke [jl] described attractors that were 
strange but not chaotic, interest in these exotic objects has been increasing 0. Several 
studies, both theoretical [Q-|J and experimental [p~0| — 15], have, over the years, elucidated the 



principal features of such attractors which appear to be generic in quasiperiodically driven 
nonlinear dynamical systems. These are geometrically strange sets (fractals) on which all 
Lyapunov exponents are either zero or negative. Further, owing to their fractal character, 
motion on SNAs is intermittent and aperiodic. 

A very important connection between classically driven systems and the Schrodinger 
equation for a particle in a quasiperiodic potential makes the correspondence between 
nonchaotic attractors and localized states []T7| — pL9f] , and this further underscores the interest 
in strange nonchaotic attractors (SNAs). These results show that the transition from an 
invariant circle to SNA in the iterative mapping is related to the transition from extended 
to localized states in the quantum system. 

Although there are many known examples of systems with SNAs, proving the existence 
of such attractors is a mathematically nontrivial task. Indeed, rigorous results exist for only 
two systems — the original system introduced by Grebogi et al., and the Harper map [|T7| . 
Keller [^] and Bezhaeva and Oseledets [2~I| have shown in the dynamical system [[I]] 



Xi+i = 2acos27r0jtanha;j (1) 
<f) i+1 = {u + (pi], (2) 

(we use the notation {y} = y mod 1) that for a sufficiently large and for u an irrational 
number, the attractor is fractal, has a singular-continuous spectrum, and is the support of 
an ergodic SRB measure p2 |. For the Harper map (E being a parameter) 



x i+ i = — [xi — E + 2a cos 27T0;] 1 (3) 
<f> i+1 = {uj + 0,}, (4) 
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a persuasive argument |L7j suggests that, again for a sufficiently large and ui irrational, the 



attractor of the dynamics is fractal, and the dynamics is nonchaotic since the Lyapunov 
exponent is nonpositive [ |T8|j . 

Virtually all known examples of systems with SNAs — and the two cases above are 
typical — have the skew-product [] form and have a quasiperiodic driving term. How neces- 
sary are these features? This question has considerable practical relevance since there are 
known cases of experimental systems where the dynamics appears to be on strange 



nonchaotic attractors, but where there is no quasiperiodic driving. Furthermore, it is also 
known that a judicious choice of forcing term can convert chaotic dynamics to nonchaotic 
dynamics [p^ -p6|. Are the resulting attractors SNAs? In particular, can SNAs be formed 
via non-periodic (and also not quasiperiodic) driving of a nonlinear system? 

In the present paper we address two issues that are of concern in the study of SNAs. 
The first deals with the ubiquity of such attractors: what are the necessary conditions for a 
dynamical system such that SNAs result? While we do not provide an exhaustive answer, 
our study shows that strange nonchaotic dynamics can be realized in a wide variety of 
systems, and may in fact be quite common. Our examples all remain, however, within the 
skew-product paradigm. The second issue relates, as discussed above, to the necessity of 
quasiperiodic forcing and we find that it is possible to construct dynamical systems wherein 
there is no explicit quasiperiodic forcing and the dynamics is on SNAs. 

These results are described in the following two sections of this paper. In Section 2, we 
extend the general arguments which establish the existence of SNAs in Eqs. ([!]) and Eq. ([3D 
so as to yield a large variety of dynamical systems wherein the motion will remain on SNAs. 
In Section 3, we further generalize the dynamics so that the forcing is no longer quasiperiodic 
but the resulting dynamics continues to be on SNAs. This is followed by a brief summary 



1 In a skew-product system, one of the variables evolves entirely independently of the other (s), as 
t> does in these examples. 
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and discussion in Section 4. 



II. STRANGE NONCHAOTIC DYNAMICS 

By now there is a host of examples of strange nonchaotic dynamics [[2J in driven systems. 
These include both flows such as the driven Duffing oscillator |27],[28] and iterative maps such 
as the driven logistic map |29]j3(| . Most studies to date have relied on numerical techniques 



to establish the presence of SNAs, by explicitly computing the Lyapunov exponents and 
by determining the fractal dimension. SNAs are "intermediate" between strange chaotic 
attractors and nonchaotic regular dynamics with regard to many dynamical and structural 
properties, some of which can be determined through study of power-spectra, correlation 
functions, and related measures (see Ref. for a recent review). The most extensively 
studied cases all have a skew-product dynamical structure, and thus these are systems in 
n+ 1 dimensions, the latter dimension pertaining to the dynamics of the driving term. 

In the present paper, we study the case of driven iterative maps in 1+1 dimension for 
convenience. Rigorous mathematical results concerning SNAs are so far available only for 
the simplest such examples. However, the arguments presented here largely carry over to the 
case of flows as well, and can be extended to higher dimensions, so this is not a restriction. 



A. Generalizing the original SNA system 

The following heuristic arguments, first put forward by Grebogi et al. [p]] for the system 
embodied in Eq. (p]) suggest that SNAs should occur for appropriate values of a. 

The mapping x — > 2a tanh x is 1-1 and contracting, taking the real line into the interval 
[—2a, 2a]. Because u is an irrational number, the dynamics in <p [cf. Eq. (0)] is ergodic in the 
unit interval. The attractor of the dynamical system Eq. ([]]) therefore must be contained 
in the strip [—2a, 2a] ® [0,1]. A point x n with the corresponding <ft n = 1/4 will map to 
(x n+ i = 0, n+ i = uj + 1/4), after which subsequent iterates will all remain on the line 
(x = 0,0), as will points with n = 3/4. This line therefore forms an invariant subspace, 
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but for large enough a, this subspace is transversally unstable. Thus, it follows that the 
attractor has a dense set of points on the line x = 0, 9 G [0, 1] (since u is irrational), but the 
entire line itself cannot be the attractor for a > 1, since the dynamics is unstable on that 
line. There will, therefore at some a, be a "blowout bifurcation" |||| transition to strange 
nonchaotic dynamics. 

Generalizing Eq. (|l|) (keeping the skew-product structure intact) as 

Xi+i = af{xi)g{(j)i) (5) 

it is clear that the same arguments will carry over so long as the following properties hold: 

(i) f(x) is 1-1 and contracting, with /(0) = 0. 

(ii) /'(0) 7^ 0, and concurrently, 

(iii) g((p) = for some = 0*. 

Then, clearly, all points x, 0* will map to 0, {u + 0*}, and the subsequent dynamics will 
be dense on the line x — 0, 0. From condition (ii) above, this can be made unstable locally 
for sufficiently large a, and from (i), since the map is contracting, the slope |/'| < 1 almost 
everywhere, the Lyapunov exponent can be made negative so as to give a SNA. 



B. A piecewise linear SNA 

A piecewise linear strange nonchaotic attractor can be simply obtained by taking, in 
Eq. ©, 

f(x) = ax \x\ < 1/a 

= sign(x) \x\ > 1/a, (6) 
g{cj>) = 0-1/2. (7) 

It is easy to verify that there is indeed a blowout bifurcation near a ~ 2.33 (Fig. la), above 
which the attractor is strange and nonchaotic: see Fig. lb. 
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C. Generalizing the Harper map 



In the Harper map, on the other hand, the argument for the existence of strange non- 
chaotic dynamics proceeds as follows [ I7f . In the strong-coupling limit, a — > oo, Eq. ([3]) 



reduces to 

x i+ i = — [2a cos27r0j~ 1 . (8) 

Clearly, for (pi =1/4 or <pi =3/4 this gives a singularity in the neighbourhood of which the 
mapping locally looks like a hyperbola. Since the (ft dynamics Eq. is ergodic in the interval 
[0,1], the image of this singularity is dense: on every 0-fiber there will be a singularity, and 
thus the dynamics is on a strange set. By continuity, even away from the strong coupling 
limit but for a large enough, this argument suggests that the dynamics can be strange. (In 
fact, for Eq. (3), there are SNAs even at a = 1.) For a suitable choice of function f(x), the 
Lyapunov exponent may turn out to be negative; in such a case, the attractor is strange 



and nonchaotic. For the Harper map ||18|| , when E = this is indeed the situation. 

Proceeding as above, one can generalize the Harper map (with the ^-dynamics un- 
changed) as 

x i+1 = [f(x l ) + 2ag((f) i )}- 1 (9) 

where / and g are now arbitrary functions, the only additional requirement being that g(<p) 
should have a zero in the interval [0,1]. For suitable functions f(x), the Lyapunov exponent 
can indeed be made zero or negative, giving therefore, SNAs. 



D. SNAs of the Fibonacci chain 



The Harper map [|T7| derives from the Harper equation |3T| which is the discrete 
Schrodinger equation for a particle in a quasiperiodic potential, 



(10) 
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ipn denoting the wave-function at lattice site n, the potential being V(n) = 2a cos 27r(na; + 
O ). The identification ip n _i/ip n = x n , connects Eq. ( |I0D and Eq. ([|). This system is known 
to have critical (or power-law localized) states for a = 1, when the classical system has 
critical SNAs jH| . It is also known that other forms of the potential V(n) support critically 
localized states B^-p5|, one example being the Fibonacci chain with 

V(n) = a 0<{nu}<uj (11) 

= -a lu < {ncu} < 1, (12) 

when all states are critical for any a. 

The classical map corresponding to this potential is 

x i+1 = -[Xi-E + Vii)}- 1 (13) 
<j) i+1 = {uj + (14) 

where V(i) is given by Eq. (|Tl"|), and there is an additional overall phase <f)Q. It is a simple 
task to compute the Lyapunov exponent for this mapping (Fig. 2a); if E is an eigenvalue, 
then the Lyapunov exponent is zero. Comparing Eq. (|9|) with Eq. (0), the conditions on 
the functions / and g are met, and therefore, the attractors of the above system are SNAs; 
see Fig. 2b for an example. 



III. SNAS WITHOUT QUASIPERIODIC DRIVING 

It is clear from the above constructions, that there are two main ingredients in achieving 
SNAs in dynamical systems such as Eq. ([5]) or Eq. ([5]). Firstly, one needs some mechanism 
for local inst ability while maintaining global stability. The main purpose of the quasiperiodic 
driving, namely the <fi dynamics which is governed by Eq. @, is to make the instabilities 
dense in 0. 

This suggests that a generalization of Eq. (||) to 

x i+1 = af(xi)g(<pi) (15) 
(f) i+1 = h((pi), (16) 

7 



or of Eq. (§) to 

= -[/(»<) + "^(^i)]" 1 (17) 
4>i+i = Hfa), (18) 

where the function h is not necessarily the rigid rotation, but is otherwise such that the 
orbit /i n (0*) is dense will still yield SNAs for sufficiently large a. 

To be mathematically more precise |36|], we need that h be a homeomorphism with an 



invariant ergodic probability measure of full support, and that the h 1 orbit of 0* be dense, 
in which circumstances, the resulting attractor can be shown to be a SNA, following the 



basic proof given by Keller [2(| for the case / = tanh, g = cos and h the irrational rigid 
rotation. The basic property that is required is that the mapping h(<f>) take the interval 
[0,27r] into some continuous subinterval (at least), and that the orbit of a typical point 
should be dense in this subinterval. 

It should be added that the function h should have only nonpositive Lyapunov exponents 
since the skew-product form for the dynamical system is being retained. There are a number 
of possible choices for h which are distinct from the quasiperiodic rotation, but which use 



related maps to generate ergodic flows. This follows from the Weyl theorem |37[ which states 
that if 7i is a polynomial function of degree r > 1 with real coefficients a , . . . , a r , at least 
one of which is an irrational number, then the map 

4> n = {7~L( n )} = { a o + aiTi + • • • + a r n r } (19) 

is ergodic, and the sequence {4> n } is uniformly distributed in the interval [0,1]. If 7i is 
nonlinear, then it can be easily verified that the sequence {0„} is not quasiperiodic, but the 
Lyapunov exponent is zero since the map preserves distance. 

An example of such a system is easy to devise. A simple choice is to take TC(n) = um 2 , u 
irrational, which gives 

x i+1 = af{xi)g{(j)i) (20) 
<f>i+i = {4>i + (2i + IV}- (21) 
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The SNA which obtains for / and g given by Eqs. (|6]) and (J7|), for a suitably large value of 
a is shown in Fig. 3a. 

Of course, it is also possible to use other mappings which generate quasiperiodic motion 
to determine the <ft dynamics. Many examples of such maps are known, as for instance the 
diffeomorphisms 

0^ {0 + fi + ea(0)} (22) 

for Q a constant, e sufficiently small and a an arbitrary analytic function. The orbits are 
everywhere dense on the interval and typically have irrational rotation number; indeed, by 
Denjoy's theorem, any orientation preserving C 2 diffeomorphism of the circle is topologically 
equivalent to a rigid rotation [fjSfl . In the specific case of the circle map with a(0) = 



sin27T0 the parameter ranges wherein the map has irrational winding number have been 
comprehensively described. 

Other possibilities for h exist: there are examples of integrable geodesic flows on compact 
manifolds which have positive topological entropy but have no positive Lyapunov exponents 
Ptjfl , or one can even take h to be a SNA map such as Eq. (Jf]) itself, since that is known to 
provide an ergodic flow pil. This yields, for instance, the system 



Xi+i = 2a0jtanhxi (23) 
<f>i+x = 2(5 cos 2ir6i tanh 0; (24) 

e t+1 = {u + e t }, (25) 

which has strange nonchaotic attractors for suitable values of a and (3, see Fig. 3b. 

IV. DISCUSSION AND SUMMARY 



Although so far all known examples of systems with strange nonchaotic dynamics have 
appeared to require quasiperiodic forcing, the present work shows that this is, in fact, not 
necessary. By suitably generalizing systems wherein SNAs are known to exist [|TT|,P~5|,[2D|,[2T 
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we have constructed dynamical systems where the motion is on strange nonchaotic attrac- 
tors and there is no quasiperiodic driving. This has particular relevance to experimental 
observations of apparently nonchaotic attractors where the dynamics does not have explicit 
quasiperiodic forcing. One specific example is of a gas discharge plasma where the light flux 
as a function of time (with current as the driving parameter) appears to lie on a SNA. This 
has been deduced by attractor reconstruction and extraction of fractal dimension and Lya- 
punov exponents jL4| . The other example |15) pertains to the effect of noise which allegedly 
reproduces the effect of many-frequency quasiperiodic driving. 

In addition, we have also shown that strange nonchaotic dynamics can be quite common. 
Our approach has been heuristic: the same deconstruction also provides a prescription 
for obtaining a large variety of dynamical systems wherein SNAs must occur. We have 
shown examples of two specific systems wherein the attractors are piecewise-linear fractals. 
Detailed analyses of such examples may prove to be simpler than the cases studied so far. 

This procedure for creating SNAs is clearly not exhaustive and there may be different 
generalizations that will also yield strange nonchaotic dynamics. In particular, we have cho- 
sen to stay within the "skew-product" class of mappings, although other possibilities |55 
may be compatible with such dynamics as well. 
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FIGURES 
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FIG. 1. The piecewise linear system with oo = (\/5 — l)/2. Plot of (a) the Lyapunov exponent 
as a function of the parameter a, showing the blowout bifurcation at a « 2.33, above which the 
Lyapunov exponent — > — oo, and (b) the strange non chaotic attractor for a=2.5 
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FIG. 2. The Fibonacci chain, where u = — l)/2. (a) Variation of the Lyapunov exponent 
with forcing parameter a. Each zero of the Lyapunov exponent corresponds to an eigenstate in 
the quantum system, and a critical SNA in the classical system, (b) A typical SNA in this system, 
for E = 0, a = 0.8326745. Note that the variable plotted on the ordinate is tanhx rather than x; 
this is merely for convenience. 
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FIG. 3. (a) SNA in the system Eq. (|20| ) where the driving is no longer quasiperiodic. The 
parameters are a = 2.5, u is the golden mean ratio, (b) SNA in the system given by Eqs. (23-25), 
where the map governing the (ft dynamics itself has SNA dynamics. The parameters a and (3 are 
both 2.5. (this figure submitted in gif format) 
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This figure "a.gif" is available in "gif" format from: 



http://arXiv.org/ps/nlin/0105011vl 



